Current noise spectrum of a single particle emitter: theory and experiment 
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The controlled and accurate emission of coherent electronic wave packets is of prime importance 
for future applications of nano-scale electronics. Here we present a theoretical and experimental 
analysis of the finite-frequency noise spectrum of a periodically driven single electron emitter. The 
electron source consists of a mesoscopic capacitor that emits single electrons and holes into a chiral 
edge state of a quantum Hall sample. We compare experimental results with two complementary 
theoretical descriptions: On one hand, the Floquet scattering theory which leads to accurate nu- 
merical results for the noise spectrum under all relevant operating conditions. On the other hand, 
a semi-classical model which enables us to develop an analytic description of the main sources of 
noise when the emitter is operated under optimal conditions. We find excellent agreement between 
experiment and theory. Importantly, the noise spectrum provides us with an accurate description 
and characterization of the mesoscopic capacitor when operated as a periodic single electron emitter. 

PACS numbers: 73.23.-b, 73.63.-b, 72.70.-fm 



I. INTRODUCTION 

The development of on-demand single electron emit- 
ters opens a promising route towards novel nano-scale 
electronics based on the coherent manipulation of only 
a single or a few electronic wave packets. The chiral 
edge channels of the integer quantum Hall effect (IQHE) , 
obtained when a strong perpendicular magnetic field is 
applied to a two-dimensional electron gas, constitute an 
ideal experimental setup to test these new concepts in the 
design of quantum electronic circuits. It is now possible 
to fabricate micron-sized electrical networks in which the 
propagation of electrons is truly one-dimensional, ballis- 
tic, and quantum coherent, thus mimicking the propaga- 
tion of photons in optical fibers. Moreover, by deposit- 
ing metallic split gates on top of the electron gas, quan- 
tum point contacts acting as tunable electron beam split- 
ters can also be implemented. Using continuous particle 
sources, single electron interferences have already been 
observed in electronic Mach-Zehnder interferometerJl^^ 
as well as two-electron interferences in similar systemgBiSl. 
However, for experiments or applications where the tim- 
ing of wave packets is important, e. g. for interference 
experiments that require particles to arrive simultane- 
ously at the scatterer, continuous sources are not useful 
as one cannot control the emission time of electrons into 
the conductor. Continuous sources then need to be re- 
placed by triggered emitterJ^Hi^ that can produce single 
particle states in a controllable and timed manner. 

A prime example of a single electron emitter is the 
periodically driven mesoscopic capacitor, consisting of 
a sub-micron cavity coupled to an edge state. The de- 
vice was first theoretically proposed by Biittiker et al^ 
who showed that the relaxation resistance is quantized in 



units of h/2e^ independently of microscopic detailfl^H^H 
Experimentally, this prediction was confirmed by Gabelli 
et al^^. In a subsequent experiment, Feve et al. showed 
that the mesoscopic capacitor, when subject to large 
periodic gate voltage modulations, can absorb and re- 
emit single electrons at gigahertz frequencies, generating 
a quantized AC currenlP. Several proposals have been 
made to coherently manipulate the single electron states 
emitted by s uch a source in order to observe two-electron 
interference^2ll2l! or electron entanglementPs'. The use of 
two-particle exchange effects has also been suggested as a 
means to visualize the single electron states generated by 
the source in a tomography protocol^^ allowing for a 
direct characterization of the interaction between a single 
electronic excitation and its environmenlpS. Coherence 
properties of the single electron states emitted by the 
source can also be analyzed by injecting particles into 
a Mach-Zehnder interferometer^. However, in order to 
facilitate these few-fermion experiments, it is necessary 
first to accurately characterize and understand in detail 
the single electron emission process. The average current 
measurements performed so far in Refs. [7] and [8] give ac- 
cess to the average behavior of the source after a large 
number of particle emissions, but are not designed to 
provide information about the elementary processes in- 
volved in a single emission event. For example, average 
current measurements cannot distinguish between the de- 
terministic emission of exactly one electron followed by 
one hole during each period of the external driving and 
a fluctuating number of emitted particles from cycle to 
cycle which still results in the same number of emitted 
charges after many periods. The statistical properties of 
the source may b e cha racterized by measuring the full 
counting statisticP^H! of electron emissions, but in this 
case as well, the short time behavior of the emitter is not 
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tions of the current noise and compare the recent mea- 
surements of the high-frequency noisd^ to two comple- 
mentary theoretical descriptions: on one hand, the full- 
fledged Floquet scattering theorji^SMl which naturally 
applies to periodically driven systems and gives highly 
accurate numerical results for the current noise, and 
on th e oth er hand, a conceptually simple semi-classical 
modeP^I^ which shows surprisingly good agreement with 
the experiment and additionally provides us with a sim- 
ple intuitive picture of the dynamical emission processes. 
Based on the excellent agreement between the measure- 
ments of the short time current-current correlations and 
the two theoretical descriptions we provide a detailed 
characterization of the single electron emitter. In partic- 
ular, we can identify parameter regimes for which nearly 
perfect and deterministic single electron-hole emission is 
achieved in every cycle. For details of the experiment we 
refer the interested reader to Ref. l34l 



FIG. 1. (color online). The mesoscopic capacitor, 
a) Schematic representation of the mesoscopic capacitor. A 
quantum point contact (with gate voltage Vg controlling the 
transparency D) couples the edge channels (red arrows) in 
the capacitor and in the electron gas. Electron/hole emission 
is triggered by the excitation drive Vexc{t) applied to the ca- 
pacitor top gate, b) Principle of single charge emission (in 
the non-adiabatic regime) using the mesoscopic capacitor. At 
® the capacitor is at equilibrium and the highest occupied 
level (HOL) is resonant with the Fermi level of the electron 
gas ep (dashed blue line). The level spacing is denoted as A. 
is electron emission phase: the HOL is promoted far above 
the Fermi energy, causing the capacitor to emit the electron 
through the quantum point contact. (J) is the hole emission 
phase: the emptied level is brought far below the Fermi en- 
ergy of the external reservoir and an electron is absorbed from 
the electron gas, i. e. a single hole is emitted. A continuous 
repetition of the sequence ®-f>(2) results in periodic emission 
of a single electron followed by a single hole. This corresponds 
to the optimal operating condition of the emitter. 



accessible. In contrast, as recently suggested by some 
of us, the waiting time distributioiP^l between individual 
charge events would provide information about the short- 
time physics, but an actual measurement of the waiting 
time distribution still remains an open and experimen- 
tally challenging task. 

We focus in this work on the short time current-current 
correlations (or frequency-dependent noise) as an exper- 
imental tool to describe and characterize the elemen- 
tary excitations generated by the mesoscopic capacitor 
when operated as a periodic single electron source. In 
two recent short papers we have separately presented 
measurement^^ and theorj^SH of the finite-frequency cur- 
rent noise of the mesoscopic capacitor. In the present 
work, we expand significantly on the theoretical calcula- 



The paper is now organized as follows: In Sec. [IT] we 
introduce the mesoscopic capacitor and describe the ba- 
sic operating principles that enable periodic emission of 
electron-hole pairs into an edge state. Sec. |III| gives an 
elaborate account of the Floquet scattering theory ap- 
plied to the mesoscopic capacitor and we discuss calcu- 
lations of the average current and the current noise both 
for a two- and a three-terminal configuration. We em- 
ploy a non-interacting model which allows us to consider 
all possible parameter ranges and operating conditions. 
Sec. |IV| compares theoretical predictions of the Floquet 
scattering theory with experimental data of the average 
current. In Sec.|V]we present a semi-classical model that 
describes the mesoscopic capacitor around the operat- 
ing conditions, where maximally one charge (electron or 
hole) is emitted in each half-cycle. The semi-classical 
model allows us to account analytically for two impor- 
tant sources of noise: when electron and hole emissions 
become rare, the current fluctuations are shot-noise like 
and the noise spectrum is white. In contrast, close to the 
optimal operating regime, where exactly one electron- 
hole pair is emitted in each cycle, the current fluctua- 
tions are dominated by the randomness of the emission 
times within a period, giving rise to phase noise with a 
Lorentzian-like noise spectrum. In Sec. |VI| we exhaus- 
tively compare measurements, Floquet scattering theory, 
and the semi-classical model, focusing both on current 
and noise in different parameter regimes including the 
shot noise and phase noise dominated limits. We discuss 
the general properties of the noise of the single electron 
emitter as well as the deviations between the Floquet 
scattering theory and the semi-classical model when the 
mesoscopic capacitor is operated away from the optimal 
conditions. Finally, in Sec. |VIl| we present our concluding 
remarks. 



II. MESOSCOPIC CAPACITOR 

The mesoscopic capacitor is depicted in Fig. [l] It con- 
sists of a submicron-sized cavity (or quantum dot) cou- 
pled to a two-dimensional electron gas through a quan- 
tum point contact (QPC) whose transparency D is con- 
trolled by the gate voltage Vg. A capacitively coupled 
metallic top-gate controls the static offset potential in 
the dot Vq, as well as the rapidly oscillating component 
T4xc(i) generated by radio- frequency excitations. A large 
perpendicular magnetic field is applied to the sample, so 
that electrons propagate along the one-dimensional chi- 
ral edge channels that form due to the IQHE. The system 
can in principle be operated at an arbitrary integer value 
of the filling factor v in the electron gas (typically v — A), 
but in any case only the outer edge channel couples to 
the quantum dot. Electrons propagating along the outer 
edge channel of the quantum dot experience a discrete 
energy spectrum with energy levels that are separated 
by a constant level spacing A, see Fig. The levels 
are broadened by the finite coupling between the quan- 
tum dot and the electron gas, determined by the QPC 
transmission D. Interaction effects within the quantum 
dot were not experimentally observed^ and are thus ne- 
glected throughout this work (the absence of Coulomb 
interactions may be due to the screening from the large 
metallic top gate as well as the presence of the inner edge 
channels in the quantum dot). 

Without the periodic driving 14xc(i) (corresponding to 
the equilibrium situation denoted as © in Fig. [T|d), the 
position of the energy levels with respect to the Fermi 
energy is determined by the constant voltage Vq applied 
to the top-gate: this top-gate voltage fixes the position 
of the highest occupied level eo (HOL) with respect to 
the Fermi energy ep at equilibrium. Adding next the 
pure AC excitation voltage Vc^c{t) causes the HOL to 
be periodically shifted up and down with respect to its 
equilibrium position. We consider the situation realized 
experimental!}'^ "^^j where a square shape excitation is ap- 
plied, causing sudden shifts of the quantum dot energy 
spectrum. The square shape excitation contains a broad 
range of Fourier components and is thus a non-adiabatic 
excitation with respect to all relevant time and energy 
scales. If the peak-to-peak amplitude of the excitation 
drive 2eVc^c is comparable to the level spacing A, the 
HOL is promoted to an energy far above the Fermi 
level in the first half-period of the drive (labeled as in 
Fig. [TJj), where the electron occupying the level is then 
emitted to the electron gas through the quantum point 
contact. In the following half-period (labeled as (2)), the 
emptied level is next brought to an energy e_ far below 
the Fermi energy, where an electron is absorbed from the 
electron gas (corresponding to the emission of a hole as 
indicated in Fig. [I]). Repeating the sequence 0o(2) at a 
drive frequency of /d ~ 1 GHz thus gives rise to periodic 
emission of a single electron followed by a single hole. 

Obviously, the discussion above depends crucially on 
the value of the static potential Vq, which fixes the equi- 
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FIG. 2. (color online). Different operating conditions for the 
mesoscopic capacitor. Only the HOL is shown; the curved 
arrow represents the AC drive (with amplitude 2e1/cxc = A) 
which periodically and abruptly shifts the HOL above and 
below the Fermi level ep (dashed blue Une). The dotted ar- 
row denotes the energy of the emitted electron measured with 
respect to the Fermi level, a) Optimal operating condition: 
the energy of the HOL without driving field Vcxc equals the 
Fermi energy of the electron gas, eo ~ ep- With the driv- 
ing field applied, the electron (hole) is then emitted at A/2 
above (below) the Fermi level ep. b) Intermediate case with 
eo = ep -I- A/4: the electron (hole) is emitted at 3A/4 (A/4) 
above (below) the Fermi level, c) Non-optimal operating con- 
dition with eo ~ ep + A/2: the electron is then emitted at ~ A 
above the Fermi level, while the hole is emitted on resonance 
with ep . In this case, emissions of spurious electron-hole pairs 
lead to unwanted electrical fiuctuations. 



librium position eo of the HOL and thus the positions 
e+ and e_ during the electron emission (0) and hole 
emission (0) phases, respectively, Indeed, as illustrated 
in Fig. [2j for certain values of eg, the HOL may be in 
resonance with the Fermi level during one of the two 
phases. Such a situation is depicted in Fig. [2];, where 
eo ~ ep + A/2 and 2eV;,xc = A, resulting in e+ « ep + A 
and e_ « ep. Thus, during the hole emission phase, the 
HOL is resonant with the Fermi level, and several charges 
can be absorbed and re-emitted during a single hole emis- 
sion phase (note that during the electron emission phase, 
the second occupied level is also resonant with the Fermi 
energy). As predicted in Ref . [351 such emissions of spu- 
rious electron-hole pairs degrade the quality of the single 
particle source and lead to unwanted electrical fluctua- 
tions. In this respect, the optimal operating conditions of 
the emitter are obtained when the HOL is alternatively 
brought far above and below the Fermi level ep during 
the emission cycle. Two such cases are shown in Fig. [2^ 
and b. However, even under these favorable conditions, a 
too large value of the QPC transmission D may broaden 
the level so much that it starts to overlap with the Fermi 
energy of the lead. The optimal operating conditions are 
therefore determined by a subtle interplay between the 
static potential Vq, the amplitude of the AC drive T^xcj 
and the transmission probability D of the QPC. 
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III. FLOQUET SCATTERING THEORY 

We now describe the Floquet scattering matrix the- 
ory used to calculate numerically the average current and 
the finite-frequency noise of the periodically driven meso- 
scopic capacitor. After a general presentation of the for- 
malism, we apply it to calculate the current and noise in 
the experimental situations considered in this work. 



A. Description of the system 

We consider the schematic setup depicted in Fig. [3] 
Electrons in the incoming edge channel can tunnel onto 
the quantum dot with the amplitude y/D = \/l — t"^, per- 
form several round-trips inside the mesoscopic capacitor, 
each taking the finite time To = l/vd, before finally tun- 
neling back out into the out-going edge state. In these ex- 
pressions, the reflection amplitude r has for convenience 
been assumed to be real and energy-independent, while I 
and Vd are the circumference of the quantum dot and the 
drift velocity, respectively. In this setup, the quantum 
dot acts as an electronic analog of a Fabry-Perot cavity. 

We first consider the simple situation, where only a 
static potential Vb is applied to the quantum dot. The 
creation w{t) and annihilation b{t) operators for an out- 
going state at time t are then related to the creation 

(t') and annihilation d(t') operators for an incoming 
state at time t' through the Fabry-Perot phase acquired 
after scattering on the quantum dot 

b{t) = j dt'U{t,t')ait'), 
U{t,t')^T5{t~t') 

OO 

9=1 

Here, eVoTo/h is the phase acquired during a single 
round-trip inside the quantum dot and e is the electron 
charge. In the Fourier domain, the creation w {e) and 
annihilation 5(e) operators for the outgoing states are re- 
lated to the operators of the incoming states a^(e) and 
a(e) at energy e ~ huj through the stationary scattering 
matrix 5(e): 

a(e) = J ^e-'^^m, 

6(e) = 5(e)a(e), (2) 

r — piro(<i-<io)/fi 
^ ' \ — j,giTo(e-eo)/n 

The density of states of the quantum dotP^ 

pie)^^S*{e)f- (3) 
Zm ae 




a(t') i^irw b(t) 
> 

FIG. 3. (color online). Schematic representation of the meso- 
scopic capacitor as a time-dependent scatterer. Electrons in 
the in-going edge channel are scattered on the quantum dot 
(red loop) subject to the time-dependent potential Vexc(i). 



consists of a series of peaks corresponding to the discrete 
energy levels of the quantum dot. The typical level spac- 
ing A = /i/to is on the order of a few Kelvins. The levels 
are broadened by the coupling to the electron gas with 
the width of the peaks given by _DA/(27r). The static 
potential eo — eVo shifts the position of the energy lev- 
els measured with respect to the Fermi energy ep which 
we thus freely can set to zero, ep = 0, throughout the 
rest of the paper. The potential shift can be written as 
a phase factor (p = e^To/h entering the expression for the 
stationary scattering matrix S{e). This allows us to de- 
scribe the position of the levels in the dot at equilibrium 
independently of the level spacing A; in particular, for 
(/) = (eo = 0), the highest occupied level is resonant 
with the Fermi energy at equilibrium, whereas for (j) = tt 
(eo = A/2), the Fermi energy lies midway in between the 
highest occupied level and the lowest unoccupied level. 

In order to induce a finite AC current in the out-going 
edge channel, one must consider a time-dependent mod- 
ulation of the quantum dot potential. In addition to 
the static potential Vb, we therefore consider a periodic 
modulation Vcxc(i) with no dc component. The period 
of the modulations is T — 1/fd — 2TT/il, which also de- 
fines the drive frequency fd- The quantum dot can now 
be viewed as a time-dependent periodic scatterer, which 
can be conveniently described using the Floquet scatter- 
ing matrix formalisrrP^. Between times t — qTo and t, an 
electron inside the quantum dot performs q round trips, 
during which it experiences the time-dependent potential 
Vexcit)- The electron then acquires an additional phasd^^ 
A(j) = J. Ve^c{t')dt'. Since the drive V^cxc(0 is peri- 

odic, we can express the acquired phase in terms of the 
Fourier coefficients c„ entering the Fourier series 



e-*t So v...(t')dt' ^ J2 c„e-™"*. (4) 

n 



The annihilation operators b{e) and a(e') then become 
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related as 

S(e) = J de'U{e,e')a{e'), 

= ^ Um{()5{e - e' + m/ifl). 

m 

The Floquet scattering theory clearly expresses how scat- 
tering occurs through the emission or absorption of a 
quantized number m of energy quanta firPH. The scat- 
tering amplitude associated with the transfer of m quanta 
is given by the Floquet scattering matrix 

t^m(e) = ^ c„c*_,_„j5(e_„), 

n 

where the notation 




FIG. 4. (color online). Experimental setups, a) Two-terminal 
geometry in which the incoming and outgoing edge channels 
(red arrows) are both connected to the same ohmic contact, 
where the current is measured, b) Three-terminal geometry 
in which the incoming and outgoing edge states (red arrows) 
are connected to ohmic contacts 1 and 2, respectively. The 
current is measured in contact 2. 



e ± mMl 



(7) 



has been introduced for convenience. We note that in 
the absence of the excitation drive, only elastic processes 
can occur and we recover U„i{e) — S{e)S„iQ in agree- 
ment with Eq. ([2]). Finally, from the unitarity of the 
time evolution operator [/, the following relations can be 
deduced: 



n 



(8) 



by the detection circuit is therefore independent of the 
parameters of the quantum dot and, in particular, of the 
QPC transmission D. The environmental noise (noise 
of the amplifier or thermal noise emitted by the detec- 
tion circuit towards the sample) is reflected on a constant 
impedance and is therefore easily subtracted in order to 
extract the noise emitted by the source as we discuss in 
detail below. Only the noise due to the capacitor is then 
measured. 



Two-terminal geometry 



B. Experimental considerations 

With the Floquet scattering matrices at hand we can 
now proceed with calculations of the average current and 
the finite-frequency noise. At this point, however, the 
experimental details of the measurement setup must be 
carefully considered. In Fig. [4] we show a two and a 
three-terminal experimental setup. In the two-terminal 
geometry. Fig. [4^, the incoming and outgoing edge chan- 
nels are connected to the same ohmic contact with a 
fixed chemical potential taken as the zero-energy refer- 
ence /i = ep = 0. The top-gate of the quantum dot 
is considered as the second terminal. The two-terminal 
setup suffices for measurements of the average current. 

However, for measurement of the current noise it is 
useful to include an additional ohmic contact as shown 
in Fig. |4]d. The additional ohmic contact is inserted be- 
tween the measurement contact and the quantum dot and 
is connected to a cold mass. In this geometry, the noise 
is measured on a constant impedance, given by the edge 
channels flowing from the measurement contact to the 
grounded contaclp]. The impedance of the sample seen 



For the two-terminal geometry in Fig. |4^, the operator 
for the current emitted by the contact is given as 

I{t) = e[&^(t)6(t) - a\t)ait)], 

de de'e'^'-''^'/''\b^{e)b{e') - a^(e)a(e') 



- / dede' e^t--^')*/?^ 
h 



X [ ^ U:,{e)Um'{e')a^{eMe'„,)-a\e)a{e') 

(9) 

Since the ohmic contact is in thermal equilibrium with 
temperature Tci and chemical potential fi — 0, we have 



{dHe)d{e'))^fie)5{e-e'), 



(10) 



where /(e) is the Fermi-Dirac distribution. The average 
current {I{t)) has the same T-periodicity as the driving 
potential Vcxdt) and can therefore be written in terms of 
its Fourier components as 



(11) 
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where 

^'^ ^ / deUKe)Ue+kie-k)[fiei) - /(e)] . (12) 
In particular, the first harmonic Iq — Ik=i is given by 

In = lY. f deUUe)Ui+i{e - hn) [f{e + mn) - /(e)] . 

(13) 

We consider next the noise emitted by the source in 
the two-terminal geometry using the expressions for the 
current operator i{t) in Eq. Since i{t) corresponds to 
a non-stationary current, the current-current correlation 
function 

cit,t')^{di{t)si{t + t')), (14) 

with 51{t) = I{t) — {I{t)), depends explicitly both on 
the absolute time t as well as the time difference t' . The 
current-current correlation function is T-periodic in the 
absolute time t, such that it can be expressed in terms of 



the Fourier components Ci{t') as 

Cit,t')^J2e'^^'Ciit'). (15) 
I 

Experimentally, the current noise spectrum is averaged 
over the absolute time t, and only the Fourier component 
Co(i') is measured. In the frequency domain we have 
Vo{uj) — J dt'e^'^* Co{t') with the mean power spectral 
density for oj > defined as 

S{u;) = 2 J dt'{5I{t)5I{t + t'))'e"^'' = 2Va{Lu). (16) 

The symbol denotes averaging over t. The noise 
for the two-terminal geometry S^ri^) contains a contri- 
bution from the cross-correlations of the current flow- 
ing from the contact, Ia{t) = ea)(t)a{t)^ and the current 
flowing into the contact, Ib{i) = {t)h{t) as well as con- 
tributions from the autocorrelations of la and /f,. The 
operators and la are related through the Floquet scat- 
tering matrix Umi^), and after some algebra we arrive 
at 



/(e)(l-/(e™ + M 



(17) 



Equations ^Vj^ and ([171 are useful for numerical calculations of the average current and the finite-frequency noise, 
respectively, for the two-terminal geometry. 



2. Three-terminal geometry 



We next consider the three-terminal geometry depicted in Fig. |4|d. In this case, electronic wave packets incident on 
the quantum dot have been emitted from contact 1, whereas the electronic waves scattering off the quantum dot are 
collected in contact 2. The total current I{t) measured in contact 2 then reads 



de de'e*('-^')*/'^ \P (e)6(e') - (e)c(e') 



(18) 



where and c are the creation and annihilation operators for edge states between contacts 1 and 2, see Fig.|4jD. Both 
contacts are in thermal equilibrium with temperature T^i and chemical potentials /ii_2 = 0, such that 



(at(e)a(e')) = (ct(e)c(e'))-/(e)<5(e-e'). 



(19) 



The average current flowing in contact 2 can now easily be calculated and we find that it is exactly equal to the 
average current obtained for the two-terminal geometry in Eqs. ( TT][l2 ). 

The current noise, in contrast, is modified compared to the two-terminal configuration. In the three-terminal 
geometry the operators for the currents flowing to and from contact 2, Ii,{t) — eb^t)b{t) and Ic{t) — ec'{t)c{t), are 
independent, such that their cross-correlation vanishes. We then find 



5'3t(w) 



2e^ 
h 



E 



de 



] C/*(e_„)C/„+™(e_„ - huj) 



/(e)(l-/(e™ + M)- 



(20) 



This expression contains a contribution from the edge channel running from contact 2 to contact 1 which is independent 
of the QPC transmission D. In order to remove this noise offset as well as any additional environmental contributions. 
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and thus only to measure the actual noise contribution from the source, we consider the difference AS3t{^) = 
S3t{^,D) — Ssri^jD = 0) between the noise at a given QPC transmission and the noise for D = 0, where the 
quantum dot is pinched off. This difference reads 



A5'3t(w) 



2e^ 
h 



E 



de 



U*{e-n)Un+rn{i-n - ^w) 



/(e) /(e + M-/(em+M 



(21) 



By construction, AS3t{<^) vanishes at zero transmis- 
sion, D = 0. Interestingly, it also vanishes at unity 
transmission, D = 1. In this case, we recover the 
noiseless flow of charges along a perfectly transmitting 
channeP^l'^. Moreover, using Eq. Q in (21), we find 
AS'3t(w = 0) = 0. The excess noise generated by the 
source has no zero-frequency component and is intrinsi- 
cally of finite-frequency natur^^H, Furthermore, it can be 
shown that AS'3t(w) = AS^ri—ojj^, such that the emis- 
sion and absorption excess noises are equal. This implies 
that S'3t(— w) — 5*37 (+w) = 2hujGK, where Gk = /h 
is the conductance of the outer edge channel flowing be- 
tween the measurement contact and the ground contact. 
This result is similar to what was found in Ref . |48] This 
implies that, as long as only excess noise AS^Ti^) con- 
sidered, no special care is required in the ordering of the 
operators entering the definition of C{t,t') in Eq. (14 1. 
One would indeed get the same result for AS'3T(w)by 
considering the inverse ordering of /(i) and I{t+t') or the 
symmetrized time-ordering. From now on, we only con- 
sider the excess noise of the source in the three-terminal 
geometry AiS'3t(w), and in order to simplify the notation 
we define S{uj) = AS^Ti^)- 



IV. AVERAGE CURRENT 
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FIG. 5. (color online). Contour plot of as functions of 
the square excitation amplitude 2el4xc/A (vertical scale) and 
the gate voltage Vg or equivalently the QPC transmission D 
(horizontal scale) . Top panel: experimental results for sample 
A from which we extract the parameters A = 4.27^, /d = 1.5 
GHz, and Td « 60 mK, see text. Bottom panel: numeri- 



cal calculations based on Eq. ( 13 1 using the extracted model 
parameters. 



The expressions for the average AC current {i{t)) 
and /n given by Eqs. ( Tlp2 ) and Eq. ([l3|), respec- 



tively, can n ow be c ompared to previous theoretical^^ and 
experimenta F * ^'^*^ works. Note that in the three latter 
references, comparisons between experimental data and 
the time dependent scattering theory were provided. In 
these models, however, the AC drive was applied to the 
ohmic contacts instead of the top gate of the quantum 
dot. This situation differs only by a simple gauge trans- 
formation when a single emitter is considered but would 
not be applicable to a system containing several emitters 
such as those investigated theoretically in Refs. and 

m 

We now consider actual experimental data: in Fig. [5^, 
we show experimental results for as functions of the 
square excitation amplitude I^xc and the QPC gate volt- 
age Vg (or correspondingly the QPC transparency D) 
taken from a sample that we label as sample A. The 
experiment was carried out in a dilution refrigerator and 
the driving frequency was fd = 1.5 GHz. The exper- 
imental results display a series of diamond-like struc- 



tures centered around 2eK:xc/A = 1. The spacing of 
the diamonds at low gate voltages Vg, where the QPC is 
nearly pinched off, allows us to extract the level spacing 
A = 4.2 K of the quantum dot. The electron tempera- 
ture Tci ss 60 mK can moreover be determined from the 
smearing of the diamond structures. Finally, the capac- 
itive coupling between the QPC gates and the quantum 
dot can be evaluated from the shift of the position of the 
levels as Vg is varied. Using these parameters we eval- 



uate numerically Eq. ( 13 1 as shown in Fig. [S] (bottom 
panel). To this end, we use a saddle-point transmission 
law for the QPCP^l. The agreement between the experi- 
mental data and numerical calculations is very good, up 
to small energy-dependent variations in the QPC trans- 
mission which were not included in the model. 

Figure [5] allows us to locate operating regimes for 
which the mesoscopic capacitor is expected to function 
optimally as a controllable single electron source. The 
white areas of the diamonds correspond to plateaus on 
which the current = 2e/d is given by the product of 
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the driving frequency and twice the elementary charge e 
(electrons and holes each contribute with an elementary 
charge, resulting in the factor of 2). In these regions, the 
device acts on average as a single electron source which 
emits in each cycle one single electron (hole) at a well 
defined energy far above (below) the Fermi level. The 
sharpness of the diamonds is then related to the accuracy 
of the expected quantization. Note that at low Vc^c these 
regions become sharp peaks in as a function of Vg, 
corresponding to the HOL being resonant with the Fermi 
energy in absence of the driving (0 — 0) (in the linear 
regime 2eVcxc hfd, keT, we also recover theoretically 
the known expression for the average AC current-^^-- ). 
In contrast, the values of V^, where two diamonds inter- 
sect, interpolate at low T4xc to zeroes in the conductance 
corresponding to the Fermi energy lying midway in be- 
tween two energy levels {4> — n). For large transmissions, 
the quantization is gradually lost and the diamonds fade 
into a linear dependence of the current on the driving 
amplitude. For low transmissions on the other hand, the 
typical escape time of an electron on the quantum dot be- 
comes much longer than the period, and charge emission 
becomes rare such that Iq is strongly suppressed. 

In the strong driving regime 2eT4xc ^ hf^, we re- 
cover theoretically well-known results for the average AC 
current'^. In particular, a second-order expansion of Eq. 
( 12 1 in the driving frequency, hil = hfd ^ -DA, enables 



us to express {I{t)) as the current response of an RC 
circuit.^ Concretely, we find 



where 



and 



7r(2fc + 1) 



-(~zl]C, + f)2i?^C2)e'(2fc+i)ot^ 

(22) 



Rq ~ 



h /c;ep^(e).F(e,Kxc) 

[/dep(e).F(e,K=xc)]' 



(23) 



(24) 



are the Vcxc-dependent capacitance Cq and resistance Rq, 
respectively, and we have defined the function 



f{e - eKxc) - /(£ + eKxc) 



(25) 



in terms of the Fermi-Dirac distribution /(e). We see 
that the capacitance Cq is given by an integral of the 
density of states p(e) over the energy window from — eV^xc 
to -t-eVexc- Under optimal operating conditions, where 
= 0, a peak in p(e) is centered around e = 0, such that 
the current becomes independent of V^xc • This is clearly 
visible on the diamonds in Fig. [5] 

In the following we restrict our considerations to exci- 
tation drives that exactly compensate the level spacing, 
i. e., 2eVc^c = A. In this case, one of the peaks in p(e) 
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FIG. 6. (color online). Time dependence of the average AC 
current {I{t)) measured in sample B for different values of the 
QPC transmission D. The level spacing and electron tem- 
perature of sample B were A — 2.5 K and T^i ~ 200 mK, 
respectively. The driving frequency was fd — 32 MHz. Ex- 
perimental data is shown with symbols, while the blue lines 
are exponential fits from which the escape time r is extracted. 
The green dashed line in panel a) corresponds to scattering 
theory calculations with 16 harmonics included. QPC trans- 
mission and extracted escape times are a) D = 0.005, r « 3.6 
ns, b) D = 0.0035, r « 5 ns, and c) D = 0.03, t 0.9 ns. 



is always fully integrated over regardless of and the 
capacitance Cq becomes independent of and Z^. In 
that case we obtain the simple result Cq — e^/A. For a 
square excitation drive voltage, the time-dependent av- 
erage current is 



-t/r 



T 1 



.-T/2t 



(26) 



for < i < T/2. Clearly, this is an exponentially de- 
caying current with a characteristic RC time given by 
the escape time t = RqCq, that is related to the QPC 
transmission D a^^ 



1 

D 



To 

D 



< 1, 



(27) 



where we recall that — h/A is the time it takes an 
electron to complete one round inside the mesoscopic ca- 
pacitor. Integrating next the current over one half-period 
of the driving we find the average transferred charge per 
half-period 



2VcxcCq X tanh(- 



^ ■) = e X tanh(— ^) . (28) 



The exponential decay of the current described by Eq. 



( 26 1 was observed experimentally in Refs. [7|and|33] Typ- 
ical experimental results are shown in Fig. [6] for a sample 
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that we label as sample jj 7 | 23 | 43 | p^j. ^^lis sample, the 
level spacing was A = 2.5 K and the electron tempera- 
ture Toi ~ 200 mK, while the experiment was carried out 
at a driving frequency of fd = 32 MHz. 

The current is well approximated by an exponential 
decay, allowing us to extract the escape time r, which is 
tunable by the QPC gate voltage Vg and thus the QPC 
transparency D. For sufficiently large QPC transmis- 
sions, the integral of the current becomes constant (see 
also data in Ref. I43p , demonstrating the quantization of 
the average transferred charge per half-period Q* — e, s o 
long as fdT ^ 1 and thus tanh(j^) ~ 1 in Eq. (28 1. 



For low transmissions, {I{t)) still exhibits an exponential 
decay with an escape time r that is comparable to the 
half-period T/2. This indicates that the RC circuit de- 
scription of the single electron emitter is still valid even 
for fdT ~ 1. The RC circuit description of the single 
electron emitter also allows us to extract t = RqCq from 
measurements of the phase of the first harmonic /n and 
Q* from both modulus and phase measurement of for 
arbitrary values of the transmission D. 

Using measurements of the modulus (see Figure[5| and 
phase of /o for sample A, the and r dependence on the 
QPC gate voltage Vg at excitation amplitude 2eVc^c — A 
are plotted in Fig. 171 As Vg is swept from large nega- 
tive voltages towards zero, the transmission D increases, 
while the escape time t decreases over two orders of mag- 
nitude. The blue line corresponds to the measurement 
data of Q* while the black dashed line expresses Q* in 



terms of r according to Eq. (28). For sufficiently short 
escape times t, Q* becomes quantized and equal to the 
electron charge e (corresponding to the quantization of 
the modulus of the first harmonic in units of 2e/ in Fig. 
[5]). Small residual oscillations around — e stemming 
from the capacitive coupling between the QPC gates and 
the cavity can be seen. When Vg is swept, the quan- 
tum dot goes periodically from optimal conditions = 
(eo = 0) to = TT (eo = A/2). In the first case, Q* is 
quantized, i. e. Q* — e, whereas in the second case, Q* 
becomes extremely sensitive to the exact value of the ex- 
citation amplitude and Q* e (slightly above e in Fig.[7| 
explaining the oscillations. The quantization of Q* can 
be checked on the inset of Fig. [Tj showing Q* as a function 
of the escape time r for conditions close to the optimal 
value of « 0. For short escape times, all data agree 
with — e within error bars (size of squares). In this 
regime, electrons and holes appear to be systematically 
emitted with the uncertainty in the emission time deter- 
mined by r. In contrast, when the escape time becomes 
comparable to half a period, the expected quantization 
is lost, Q* < e, reflecting that single charges are not de- 
terministically and periodically emitted from the quan- 
tum dot. Furthermore, Q*'/e (squares) is well described 
by the expected tanh(l/4/dT) dependence (continuous 
black line). Finally, we note that the Floquet scattering 
matrix theory predicts rapid oscillations in the currenlpSl 
on time scales on the order of Tq = h/ A which is much 
smaller than the periods considered here. The observa- 
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FIG. 7. (color online). Measured average transferred charge 
per half-period Q* (blue line) and escape time r (red line) for 
sample A as functions of Vg. Parameters are A = 4.2 K, Td ~ 
60 mK, and fd = 1.5 GHz. The dashed line corresponds to Q' 
determined from the measured r using Eq. (281. The dotted 
circle corresponds to the point at which r = T/2. Inset; 
average transferred charge for optimal operating conditions 
(0 ~ 0) as a function of the escape time r in units of T/2. 



tion of such oscillations, however, is not yet within ex- 
perimental reach, as the time scale Tq ~ 10 ps is still well 
below the measurement resolution of our experiment. 

As we have seen in this section, the measurements of 
the average AC current suggest that the mesoscopic ca- 
pacitor in certain parameter regimes acts as a control- 
lable single electron source. Average measurements only, 
however, do not reveal possible fluctuations of the quan- 
tized current and to this end we need to consider the 
noise properties of the mesoscopic capacitor. As already 
discussed, the noise spectrum can be calculated numeri- 
cally using the Floquet scattering matrix theory. In order 
to understand in detail the noise properties of the meso- 
scopic capacitor we flrst analyze a simple semi-classical 
model of the charge transport. The description of the 
mesoscopic capacitor in terms of an AC driven RC cir- 
cuit forms the basis of the semi-classical model described 
in the following section. 



V. SEMI-CLASSICAL MODEL 

As mentioned above, the noise of the mesoscopic ca- 
pacitor can be calculated numerically using the Floquet 
scattering theory. However, except for certain limiting 
cases, it is generally difficult to obtain analytic results 
from which one could hope to develop a deeper under- 
standing of the noise properties. In the adiabatic regime, 
where the period of the driving is much longer than 
any other time scale, the noise spectrum has been found 
analyticallj'2SH38| contrast, in the situation that we 
consider here, where the driving potential is highly non- 
adiabatic, a perturbative expansion in the driving fre- 
quency is not possible. Nevertheless, both the experi- 
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mental and the numerical results obtained from the Flo- 
quet scattering approach suggest that it may be possible 
to calculate analytically the noise spectrum using a con- 
ceptually simple semi-classical model. The model that 
we now describe was first suggested by Mahe et al^ and 
later investigated theoretically by Albert et alW^. In the 
following we discuss the semi-classical model from which 
we derive the finite-frequency noise spectrum and provide 
a thorough comparison between analytic, numerical, and 
experimental results. 

The semi-classical model can be explained by consider- 
ing again Fig. [TJd. The model assumes that the quantum 
dot can emit at most one electron and one hole per period 
and time is discretized in units of Tq; the time is takes 
an electron to complete one round inside the mesoscopic 
capacitor. In the emission phase @ the probability in 
each time step for an electron to escape is equal to the 
transparency of the QPC, namely D. Additionally, since 
the amplitude of the periodic driving is on the order of 
the level spacing A, higher-lying states can safely be ne- 
glected and maximally one electron can escape the meso- 
scopic capacitor as re-filling is not possible. Similarly, in 
the absorption phase (2) the probability of emitting a hole 
in each time step is D. This semi-classical model can be 
theoretically formulated as a master equation in discrete 
time for the probability of the mesoscopic capacitor to 
be occupied by an electron. Setting the electron charge 
e = 1 in the following, this probability is equal to the 
average (additional) charge of the mesoscopic capacitor 
(Q), where Q = 0,1. 

The master equation determines the evolution o f the 
average charge after one time step and takes the fornPHMl 

r {I ~ D){Q{tk.i)) ® 

{Q{tk+i,i)) = { , (29) 

[ D[l-{Q{tk,i))] + {Q{tk,e)} ® 

where we have used that 1 — (Q) is the probability for 
the mesoscopic capacitor to be empty and t = tk/ de- 
notes time at the fc'th time step during the rth period. 
The emission (absorption) phase ® ((2)) corresponds to 
k = 1,2...,K {K + 1,K + 2,..., 2K), where K is the 
number of time steps in the absorption and emission 
phases, respectively, each of duration T/2. 

Experimentally, the noise measurement frequency lu 
was roughly equal to the driving frequency ft = 27r/d and 
both were much smaller than the inverse round trip time 
To, i. e. 27r/a; ~ 27T/ft ~ 60rj221. This allows us eventu- 
ally to consider the continuous-time limit of the model, 
where the step size becomes irrelevant and drops out 
of the problem. Interestingly, the physics of the system 
is then governed by the single dimensionless ratio T/r of 
the period T over the escape (or correlation) time 

(30) 



Hi/{i-D)Y 



In the limit D ^ 1, we recover the expression in Eq. 
(27). At the end of this section, we discuss the physical 



The master equation can be understood by consider- 
ing the average current as it was calculated using Floquet 
scattering theory in Ref. |3H1 For the square-shaped driv- 
ing considered in this work, the current consists of one 
step-like term with time step Tq contained in an expo- 
nential envelope function and one oscillatory part with 
period Tq. The latter corresponds to the rapid oscilla- 
tions of the current mentioned in the previous section. 
These oscillations are due to quantum interferences be- 
tween different orbits in the mesoscopic capacitor and 
they vanish at high temperatures. Still, at arbitrary tem- 
perature only the first step-like term survives after inte- 
gration over the time step Tq, which leads to the mas- 
ter equation. The semi-classical description reflects that 
the oscillations due to quantum interferences are irrele- 
vant for the average current on a time scale that is larger 
than To- 

At this point we have not provided a detailed deriva- 
tion of the model which would require us to compare 
not only average current but also noise and higher-order 
correlations with the full quantum theory. This still re- 
mains a challenging and open task and for now we simply 
rely on the excellent agreement with experimental data 
as we demonstrate in the following. Obviously, the semi- 
classical description cannot be correct under all operating 
conditions and already now we can anticipate situations 
where the model will differ from the full Floquet scatter- 
ing theory: for example, as the noise measurement fre- 
quency approaches the internal frequency of the quantum 
dot, namely the frequency h/A = I/tq associated with 
the level spacing A, we expect that the semi-classical 
model will no longer be valid. The model also neglects 
the possibility of emitting two electrons within the same 
period and will therefore not apply to situations, where 
a level of the quantum dot is in resonance with the Fermi 
level of the 2D electron gas, i. e. for <j) — tt. We make 
a detailed comparison between the semi-classical model 
and numerics in the following section. 

Before turning to calculations of the noise spectrum, 
we consider the average charge (Q) in the mesoscopic 
capacitor. In Ref. .35, the average charge was used to 
obtain the average current flowing out of the mesoscopic 
capacitor (/(i)) = -(Q(i)) ^ [{Q{t)) - {Q{t + To))]/to. 
Solving Eq. (29) for the average charge (Q) we readily 
find 



{Q{tk,i)) = 



-(tfc,«-<!T)/T 



® 



(31) 



Here, we have defined e 



-T/2t 



ae - 1/(1 + e)+ 9e 



meaning of the correlation time. 



2£ 

and f3( — 1/(1 + e) — 9e'^^^^ with 9 depending on the 
initial conditions at the time when the periodic driving 
is turned on. The correlation time r determines the time 
scale over which the system loses memory about the ini- 
tial conditions encoded in 9 and (Q) becomes periodic in 
time. We notice the close relation of the above expres- 
sion with the average current of an RC circuit. Impor- 
tantly, the model reproduces the expressions in Eqs. (26) 
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FIG. 8. (color online). Current correlators within the semi- 
classical description, a) Shot noise limit: for small escape 
probabilities D, charges are not systematically emitted within 
each half-period: the emission probability per half-period 
P = Q'/e is small, and the current fluctuations correlator 
is essentially given by the Dirac peak in Cq'^^I'). The cor- 
relators C^"\t'), cP^^\t') and Co{t') for long time t' are 
shown in the lower panel (with black, red, and blue lines, 
respectively): Co(t') is negative for t' < 3T, reflecting the 
anti-bunching of emitted charges, b) Phase noise limit: for 
sufficiently large escape probabilities (here, D — 0.2), charges 
are systematically emitted {P ~ 1), and Cq'^''^^ {t') consists of 
a peak with a finite width given by the escape time r. In this 
case, Co(t') is negative in a smaller range of times t' < T/2. 



and ( 28 ) for the time-dependent average current and the 
average transferred charge per half-period, respectively, 
obtained from the Floquet scattering theory. 



VI. NOISE SPECTRUM 

We are now ready to discuss the noise properties of the 
mesoscopic capacitor. As we will see, the finite-frequency 
noise spectrum allows us to characterize the mesoscopic 
capacitor as a single electron source as well as to deter- 
mine the optimal operating conditions of the device. Be- 
fore presenting any detailed calculations we discuss the 
two primary sources of noise. 

A. Sources of noise 

In the semi-classical model, the mesoscopic capacitor 
emits at most one electron and one hole per period. Still, 
depending on the ratio between the escape time r and the 
period T, the source may fail to emit. We quantify the 
emission probability per half-period by the ratio 



having recalled that Q* is the average transferred charge 
per cycle, see Eq. (28 1. We refer to the noise associated 



with such cycle missing events, where the mesoscopic ca- 
pacitor fails to emit, as shot noise. However, even when 
the emission probability is close to unity and the meso- 
scopic capacitor emits an electron and a hole in almost 
every cycle, there are still fluctuations in the current as- 
sociated with the random emission times within a period. 
This source of noise is referred to as phase noise. The two 
uppermost panels of Fig. |8] illustrate the main sources 
of noise by showing typical realizations of the current, 
where emissions of electrons (holes) are shown with filled 
(empty) circles, on top of the average current. The up- 
permost panel of Fig. [8^ illustrates shot noise, while the 
uppermost panel of Fig. [Sja corresponds to phase noise. 

From the definition of the current-current correlation 
function in Eq. (|14|) we can immediately write 



(33) 



c{t,t') = {mi{t+t'))-{m){i{t+t')), 



where I(t) is no longer a quantum mechanical operator, 
since we are considering a semi-classical description. As 
already mentioned, the system is not translational invari- 
ant in time due to the periodic gate voltage modulations. 
The correlation function therefore does not only depend 
on the time difference t' , but also on the absolute time 
t. Experimentally, the correlation function is averaged 
over the absolute time t and the time-average correlation 
function is then 



(34) 



where C^"\t') = {I(t)I{t + f))^ and C'o"''it') = 

{I{t)){I{t + t')) . The time-average correlation function 
is exactly the Fourier component Co{t') entering Eq. ( 15 ) 
and the corresponding noise spectrum is consequently 
given by Eq. ^16^ . 

Figure |8] shows numerical calculations of the correla- 
tion functions in the two limiting cases. The results were 
obtained in numerical simulations of the stochastic pro- 
cess defined by the semi-classical model. The left pan- 
els correspond to the shot noise regime, where the es- 
cape time T is much larger than the period T, whereas 
the right panels show results for the phase noise regime. 
Depending on the escape probability D [or equivalently 
the rati o between the escape time and the period, see 
Eq. ffl], the contributions from C^"\t') and C^^''^^\f) 



P = Q'/e, 



(32) 



vary. The two central panels of Fig. [8] show the correla- 
tion functions Cq^^^I') and c'^^^^^\t') for short times 
t' <€i T. In both regimes, Cg^^^ (t') contains a Dirac peak 
at t' = 0. Indeed, since maximally one charge is emit- 
ted per half-period, the short-time correlations vanish. 
In this respect, the Dirac peak is the hallmark of single 
particle emission. Correlations are only recovered when 
t' is close to a multiple of the half-period, as seen in the 
lower panels of Fig. [§) 

The height of the Dirac peak at i' = is propor- 
tional to the average transferred charge per half-period 
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0*: C^"\t' = 0) counts the average number of peaks 
and dips in the instantaneous current I(t) corresponding 
to emitted electrons and holes. The correlation function 
cjf'^^^'^ it') is given by the autocorrelation of the exponen- 
tially decaying average current. At short times, it there- 
fore has a peak &i t' — Q with a finite width given by the 
escape time t. At times comparable to multiples of half a 
period, correlations are again recovered, which compen- 
sate the long-times correlations in c'^^^\t'). Thus, for 
long times t' ^ T, the correlation function C(i') van- 
ishes since charges emitted by the source are no longer 
correlated. The timescale on which C{t') decays to zero 
depends on the transmission D and is given by the escape 
(or correlation) time r. 



1. Shot noise 

For small escape probabilities 13 <C 1, the escape time 
T becomes much larger than the period T. The peak 



at = in C, 



(-f>(^) 



{t') correspondingly becomes much 



smaller than the Dirac peak in CQ^\t')^ see Fig.js^i. The 
current correlator C{t') is thus given by a Dirac peak 
at = combined with negative values up to at finite 
times on the order of the correlation time r. The noise 
power spectral density is then constant, except at zero 
frequency, where it vanishes because the integrals over 
CQ^^\t') and CQ^^^^\t') cancel each other. In this case, 
the source randomly emits charges and the charge fluc- 
tuations are similar to shot noise. 

The negative values of C(t') at finite times reflect the 
anti-bunching behavior of emitted charges: at low trans- 
mission probabilities D, the anti-bunching extends over 
a large range of times, t' ~ r, see Fig. [8^ (lower panel). 
Even in the shot noise regime, a hole must be emitted 
after the emission of an electron before a second electron 
can be emitted. Now, approximating Cg^^^^^(t') ~ 0, 

see Fig, [sk, and writing dJ'\t') = {2fdQ* /e)d{t'), we 
find* 



e 

5'shot('^) =4eVd X F = — , 

r 



(35) 



where P the emission probability per half-period, Eq. 
(32). Interestingly, this expression is identical to the 
usual shot noise formula S = 2e/p, where the usual 
charge current has to be replaced by the particle current, 
Ip, given by the sum of the average number of electrons 
emitted in the first half period and holes in the second 
one, times the product of the electric charge with the 
drive frequency: /„ = 2efd x P. 



2. Phase noise 



the emission probability P is close to unity. The time- 
dependent average current then consists of well-defined 
exponential decays with a decay time given by the escape 
time T < T/2: {I{t)) = ±(e/r)e-*/^. Here the different 
signs correspond to the emission of electrons or holes. In 
this case, we find a simple expression for the correlation 
function 



it') = ^e-l*'!/^ 



and the noise is then given by 



'S'phasc(w) = 4e^/d 



I34I35I 



1 +Cj2t2 



(36) 



(37) 



Even if charges are systematically emitted each period, 
we find a finite value of the noise which depends only on 
the escape time r. This noise is due to the uncertainty 
in the emission time of charges within a period, and is 
thus referred to as phase noise. The phase noise is an 
intrinsically high-frequency noise and it is the signature 
of single charge emission: when the source periodically 
emits single charges, the noise reduces to the value of the 
phase noise determined only by the temporal extension 
T of the emitted wave packets. 



B. Analytic expression 

We now present an analytic calculation of the noise 
power spectrum which covers both limiting cases as well 
as the intermediate regim^^. To this end, it is useful 
to consider the charge correlation function CQ{t,t') = 
{Q{t)Q{t')) - {Q{t)){Q{t')) together with the relation 
Vo{ijj) — uj'^Vq{uj). Here the definition of 'Pq{uj) is simi- 
lar to that of 7^0 (i^) in Eq. ( jlG] ), but with the current / re- 
placed by the charge Q. The charge correlation function 
can be evaluated following the schematic illustration in 
Fig.[9) We first note that {Q(t)Q{t+t')) is the joint prob- 
ability for the capacitor to be charged with one electron 
both at time t and at time t-\-t' . Using conditional proba- 
bilities we then write {Q{t)Q{t + t')) = {Q{t)){Q{t + t')), 
where {Q{t + t')) is the probability that the capacitor 
is charged with one electron at time t -\- t' given that 
it is charged at time t. For t' > 0, the conditional 
probability {Q{t + t')) can be found by propagating for- 
ward in time the condition {Q{t)) = 1 using the master 
equation in Eq. (29), see also Fig. [oj Similar reasoning 
applies to the case t' < 0. Finally, integrating over t, 
the time-averaged charge correlation function becomes 

{6Q{t)6Q{t + t'))* = J e-l*'l/^tanh {^), and we imme- 
diately obtain the noise power spectrum a^^ 



S{uj) = Ae^fa tanh 



(38) 



In the phase noise regime, the escape probability D is 
so high that charges are emitted in nearly every cycle and 



Interestingly, the noise power spectrum is given by 
the average charge emitted during the emission phase 
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<Q(t)> A 



<Q(t)> <Q(t)>^l 



<Q(t+t')> 







r 



r 



— 



FIG. 9. (color online). Schematics of periodic driving and (conditional) charge occupations. The mean occupation of the 
mesoscopic capacitor {Q{t)) is shown in red. The mean occupation is equal to the probability of the mesoscopic capacitor to 
be occupied with one (additional) electron. The blue curve shows the conditional mean occupation {Q{t + t')). This is the 
probability that the mesoscopic capacitor is charged with an (additional) electron at time t + t' , given that it was charged at time 
t, such that {Q{t)} = 1. Correlations decay on a time scale set by the correlation time r, implying that {Q{t + t')} ~ {Q{t + t')) 
for t' > r. 



tanh(l/4/(;r) (the factor of 2 accounts for the additional 
contribution from the average charge absorbed in the ab- 
sorption phase) multiplied by a Lorentzian-like frequency 
dependence, which accounts for the exponential decay of 
correlations in the time domain with time constant t. 
Finally, the factor w^r^ reflects that the noise spectrum 
becomes flat in the high-frequency limit, while the zero- 
frequency limit S{0) = shows that charge on average 
does not accumulate on the capacitor once (Q) has be- 
come periodic in time. From Eq. ( 38 ) it is straightforward 



to recover the limiting cases given by Eqs. (35) and (37 1, 



but the analytic result above also accounts for the inter- 
mediate regime where the escape time r is comparable 
to the period T. 



C. Detailed comparison 

We are now ready to carry out a careful comparison 
of experimental and theoretical results. In Fig. [TO] we 
first compare results for the noise spectrum obtained 
from the full Floquet scattering theory and the semi- 
classical model. We focus here on the experimentally 
relevant regime with w « $7 <C A/S together with the 
optimal operating conditions of the source, cj) = 0, and 
consider three different values of the transmission prob- 
ability D. The figure shows that the two complementary 
approaches yield results that are in excellent agreement. 
We observe that our numerical and analytic calculations 
agree well both in the shot noise limit with D — 0.01 
(t « 3.5 T/2) and in the phase noise limit with D = 0.2 
(t « 0.16 T/2), as well as in the intermediate cross- 
over regime. For ujt ^ 1, the phase noise S'phasc sat- 
urates to ie^fd, such that S{uj) ~ Ae'^fdP as seen both 
for D = 0.06 (r w 0.57 T/2, P = 0.71) and D = 0.01 
(P = 0.14). 

The two theoretical approaches can also be compared 
with experimental data, obtained with sample A using a 



CM 

(D 2 



semi-classical model : scattering theory, (p=0 : 



-D=0,01 
D=0,06 
-D=0,2 

■ phase noise limit 

■ shot noise limit 



H □ D=0,01 
D=0,06 
D=0,2 




FIG. 10. (color online). Finite- frequency noise S{u)) close to 
the driving frequency, a; ~ f2. Symbols; scattering theory 
((/) = 0); continuous lines: semi-classical model; dashed lines: 
limits given by Eqs. (351 and (37 1. Parameters are A = 4:K 
and fd = 1.5 GHz. 



high sensitivity microwave noise measurement setup im- 
plemented in a dilution refrigeratoi'^. The noise was 
measured in a bandwidth of 1.2 — 1.8 GHz centered 
around the driving frequency fd — 1.5 GHz. Figure [TT] 
shows the dependence of the noise ^(a; = ft) on the es- 
cape time r. The numerical and analytic results (over- 
lapping continuous lines) are in excellent agreement with 
the experimental data (circles) obtained at (/) « 0. In par- 
ticular, the experimental results are captured by the an- 
alytical expression in Eq. ( 38 ) over more than two orders 



of magnitude of the escape time t, going from the phase 
noise limit (blue dashed line), where the source exactly 
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T(T/2) 1 



FIG. 11. (color online). Noise power spectrum S{ijj = Q) as 
a function of the escape time r. Left axis: measured average 
transferred charge per half-period Q* (black squares) together 
with Eq. ( |28[ l (continuous black line). Right axis: measured 
noise spectrum (red circles, « 0, sample A) together with 
Eq. ( |38[ ) (re d lin e). Dash ed lines correspond to the two noise 
limits, Eqs. (35 1 and (37 1, illustrated schematically in the top 
panels. 



scattering theory: / ^ 
□ A=2K //'^\\ 
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shot noise limit 

phase noise limit 
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FIG. 12. (color online). Universality of the noise spectrum. 
Noise spectrum based on the semiclassical model, Eq. 1 39 1 
(blue line), and full numerical Floquet scattering theory for 
different values of the level spacing (symbols) as a function 
of the escape time r with Td = 100 mK. The dashed lines 
correspond to the two limits given by Eqs. (1351) and (|37| 



1. Universahty 



emits a single electron and a single hole in each cycle, 
to the shot noise limit (black dashed line) where particle 
emissions are rare and shot noise like^^. As discussed in 



Sec. IV the average emission probability P (squares) is 
also well captured by the expected tanh.1/4 fdT depen- 
dence (continuous black line). The results shown in Fig. 
[TT] demonstrate that the mesoscopic capacitor indeed be- 
haves as an controllable on-demand single electron emit- 
ter when operated under optimal conditions. 



Figure [TO] clearly illustrates that the Floquet scatter- 
ing theory and the semi-classical model agree well under 
optimal operating conditions. Moreover, by plotting the 
noise power S{uj — Cl) in units of the driving frequency 
e^/rf as a function of the escape time r in units of the 
half-period T/2 one finds strong indications of a simple 
universal behavior of the noise power which is indepen- 
dent of the specific parameters of the system. Indeed, by 
rewriting Eq. (38) in term of these normalized units, we 
find 



D. Universality and deviations 



^2 



While the semi-classical model is only valid for a re- 
stricted range of parameters, the Floquet scattering the- 
ory in contrast allows us to explore the full set of experi- 
mentally relevant operating conditions, including changes 
of temperature, level spacing, and measurement fre- 
quency uj/2'K. In the following, we first discuss a particu- 
lar universal property of the noise under optimal operat- 
ing conditions as described by the semi-classical model. 
Secondly, we discuss possible discrepancies between the 
semi-classical model and the full numerical calculations 
based on the Floquet scattering theory. As we will see 
below, the two descriptions start to differ once the sys- 
tem is not operated under optimal conditions, or when 
the noise is measured at very large frequencies huj ^ A. 



from which this universality is evident. In particular, 
neither the level spacing A nor the temperature enter 
this expression. The universality of the noise spectrum 
is verified by our full numerical Floquet scattering theory 
calculations with different values of A as shown in Fig. 
[12] (the dependence of r on A is taken into account). We 
also checked performed calculations with varying temper- 
atures (not shown) and found good agreement with the 
expression above. The universal behavior can be under- 
stood by noting that under optimal operating conditions, 
the noise arises from elementary charge transfer processes 
which only depend on the parameter r. As long as the 
charges are emitted sufficiently far above or below the 
Fermi level, these processes do not depend of the energy 
at which the charges are emitted. 
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FIG. 13. (color online). Noise spectrum as a function of the 
QPC gate voltage Vg. Circles: experimental data (sample 
A); continuous line: Floquet scattering theory (A = 4 K, 
Toi — 100 mK, 3 odd harmonics in the excitation drive). The 
panels illustrate charge emission in the minima {(j> = 0) and 
the maxima (0 — n) of the noise oscillations. 



2. Deviations 

Deviations from the universal behavior start to appear 
as the QPC transmission approaches unity, D ~ 1, and 
the escape time r becomes comparable to the inverse 
level spacing Tq = h/A. For short escape times, Eq. (39) 
predicts that the noise would vanish. However, the semi- 
classical description is expected to break down as the 
relevant time scales of the problem approach Tq. Small 



deviations between the semi-classical model and full nu- 
merics are already visible in Fig. 12 for r < 0.1 T/2. 



The semi-classical model is also not expected to be valid 
when one of the levels in the quantum dot is brought into 
resonance with the Fermi energy during the emission cy- 
cle. In this case, the total charge on the quantum dot 
is no longer quantized and the quantization of the first 
harmonics of the average AC current is lost, see Fig. [5] 
Under these conditions, the noise is expected to depend 
strongly on various parameters such as temperature, the 
shape of the excitation drive, and the static potential in 
the quantum dot. 

The expected parameter dependencies are clearly vis- 
ible in the noise measurements. In Fig. [13] we show the 
measured noise spectrum as a function of the QPC gate 
voltage Vg (red circles), which simultaneously controls 
the QPC transmission D and the levels in the quantum 
dot via a capacitive coupling, see also Fig.[5j Oscillations 
in the noise are observed for —0.345 V < < —0.330 V 
(0.8 ^ f ^ 1). The maxima of the oscillations cor- 
respond to the case = tt, where one of the levels is 
brought into resonance with the Fermi energy, while the 
minima occur at the optimal operating conditions = 0, 
see insets in Fig. [T3j The semi-classical model cannot 
account for these oscillations-''. Instead, we have used 
the Floquet scattering theory and numerically evaluated 
Eq. ( 21 1 as a function of Vg. The dependence of the QPC 
transmission on the gate voltage D{Vg) was extracted us- 
ing Eq. (27) in combination with the escape time r as a 




FIG. 14. (color online). Noise spectrum as a function of the 
measurement frequency lj. Blue dotted line: semi-classical 
model. Since time is discretized in units of To = h/A, S{uj) 
can only be calculated up to huj = A/2. Black line: Floquet 
scattering theory with (f) = 0. Red dashed line: Floquet scat- 
tering model with = tt. The electronic temperature was 
set to Tel = 100 mK. The blue circle denotes the noise at the 
driving frequency, uj — Q,. Parameters are D = 0.2 {P ~ 1), 
A = 4 K, and fa = 1.5 GHz. 



function of Vg, Fig. ItI The AC drive Kxc(i) used in the 
calculations consistea of three odd harmonics, accounting 
for the finite bandwidth of the microwave pulse generator 
used in the experiments. 

are 



The numerical results (blue line) shown in Fig. 13 



in good agreement with the experimental data. The max- 
ima in the noise can be understood by noting that a 
square-shaped AC drive with a finite number of harmon- 
ics contains fast oscillations, or ripples, which affect the 
energy resolution of the emitted charges. In the case 
(j) = TT, the level brought into resonance with the Fermi 
energy and then oscillates rapidly with respect to the 
Fermi energy. This "shaking" of a resonant level causes 
additional charge transfers (which also depend on the 
QPC transmission), leading to an increase in the noise. 
Due to such spurious emissions of electron-hole pairs, the 
source does not behave as a perfect single particle emitter 
under these operating conditions. 

Deviations between the Floquet scattering theory and 
the semi-classical model also occur when the measure- 
ment frequency becomes comparable to the level spac- 
ing u! ~ A/h. Indeed, since the semi-classical model 
describes the dynamics of charge emissions with a dis- 
crete time step Tq — h/A, it cannot account for dynam- 
ics on time scales shorter than r, or large measurement 
frequencies w ~ A/h. The Floquet scattering model, on 
the other hand, predicts fast oscillations on time scales 
comparable to Tq in the average AC current '^^ . We there- 
fore expect strong deviations between the semi-classical 
model and the Floquet scattering theory in the noise 
spectrum at high frequencies. 

Numerical calculation of the noise spectrum at high 
measurement frequencies are shown in Fig. [14] obtained 
from the Floquet scattering theory and the semi-classical 



16 




o o o 



<n>/N 



<An>/<n> . 

0.2 0.4 0.6 0.8 
D 



FIG. 15. (color online). Mean and variance of the counting 
statistics of emitted electrons P{n, N) after a large number 
of periods A'^ as functions of the QPC transparency D. Full 
lines: Results based on the semi-classical model. Symbols: 
Numerical data obtained from the Floquet scattering theory. 
Numerical results (symbols) were adapted from Ref. 1271 



model. We see that the two complementary approaches 
agree well at small frequencies for = 0. However, 
while the noise in the semi-classical model saturates to 
S{uj V,) = 4e^/d at high frequencies, the Floquet scat- 
tering theory, in contrast, is cut off for huj w A/2, where 
it drops to zero. Indeed, the electrons (holes) are emit- 
ted at an energy A/2 above (below) the Fermi energy, 
corresponding to emission of radiation (or photons) at 
frequencies below A/2. For (p = tt, charges can be emit- 
ted at energies up to A above or below the Fermi energy, 
and the cutoff frequency is then equal to A as seen in the 
figure. We note that the excess noise shown here indeed 
is symmetric in uj as expected. 



E. Full counting statistics 

We round off this section by discussing the counting 
statistics of emitted electrons. Under optimal operating 
conditions, the semi-classical model fully accounts for the 
charge dynamics of the emitter at low frequencies and it 
allows for tractable calculations of the counting statistics 
P(n, N) of the number of emitted electrons n during a 
large number of periods. In principle, we can calcula te all 
moments (or cumulants) of the distributiorP^'^^Ill]^ l^ut 
we focus here on the first and second cumulant, the mean 
(n) and the variance (An^) with An = n — (n). In Fig. 
15 we show a comparison between the calculations of the 
first two cumulants based on the semi-classical modeP^ 
and Floquet scattering theorjEH. We observe an excel- 
lent agreement for small transparencies, but eventually 
deviations appear for D > 0.3. These discrepancies ap- 
pear as the broadening of energy levels become so large 
that the effect of spurious emissions of (several) electron- 
hole pairs during a period becomes non-negligible. In 



this case, the mean number of emitted electrons during 
a period can exceed one. The counting statistics is im- 
portant for characterizing the accuracy of the mesoscopic 
capacitor as a single electron emitter^. 



VII. CONCLUSIONS 

We have investigated experimentally and theoretically 
the finite-frequency noise spectrum of the mesoscopic ca- 
pacitor when operated as a periodic single electron emit- 
ter. We have compared experimental data with two com- 
plementary theoretical descriptions. On one hand, we 
discussed the Floquet scattering theory which allows us 
to accurately describe the system over the full range of 
experimentally relevant parameters, in particular the en- 
ergy ranges in which charges are emitted. On the other 
hand, we considered a semi-classical model which, despite 
its simplicity, is able to account for the charge dynamics 
of the emitter when operated under the optimal operating 
conditions. This model allowed us to develop an analytic 
understanding of the measured noise spectrum and the 
numerical results obtained using the Floquet scattering 
theory. 

Depending on the escape time of electrons from the 
mesoscopic capacitor, two distinct noise regimes could be 
identified. When the escape time is much smaller than 
the period of the drive, the mesoscopic capacitor emits 
electrons and holes in an almost fully periodic manner 
with unit probability and the main source of noise is due 
to the uncertainty in the emission time within a period. 
This type of noise, referred to as phase noise, can be 
clearly identified both in the theoretical and experimen- 
tal results. The phase noise provides a fundamental lower 
limit on the noise, arising from the random jitter be- 
tween the triggering of emission and the actual emission 
time. In the phase noise regime, we obtained excellent 
agreement between our experimental data, the Floquet 
scattering theory, and the semi-classical model. In the 
other extreme, when the escape time of electrons is much 
larger than the period of the drive, electron emission be- 
comes rare and the fluctuations are shot noise-like with 
a white spectrum related to the average (particle) cur- 
rent Ip = 2efd X P through the usual Schottky formula 

S'Schottky = 2elp. 

As the mesoscopic capacitor is tuned away from the op- 
timal operating conditions and charges are emitted close 
to the Fermi energy, a significant increase in the noise is 
observed due to additional charge fluctuations generated 
by the source. These spurious emission processes are not 
accounted for by the semi-classical model, in which max- 
imally one electron-hole pair can be emitted during each 
cycle. In contrast, these additional fluctuations are fully 
accounted for by the Floquet scattering theory. The abil- 
ity to accurately investigate, model, and characterize the 
single electron emission process will prove useful in fu- 
ture few electron experiments involving the mesoscopic 
capacitor as a controllable single electron source. 
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